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Abstract 

We derive a general quantum formula giving the mean-square displacement 
of a diffusing particle as a function of time. Near K we find a universal 
logarithmic behavior (valid for times longer than the relaxation time), and 
deviations from classical behavior can also be significant at larger values of 
time and temperature. Our derivation depends neither on the specific compo- 
sition of the heat bath nor on the strength of the coupling between the bath 
and the particle. An experimental regime of microseconds and microdegrees 
Kelvin would elicit the pure logarithmic diffusion. 
The so-called fluctuation-dissipation theorem-which relates the thermal fluctuations of a 
variable x to the response of that variable to a weak external force-is usually described as 
generalizing the Smoluchowski- Einstein relation for Brownian motion, D = kT{x; but it is 
not easy to find in the literature any explicit derivation of this relation as a direct corol- 
lary of the theorem. In this paper we will provide such a derivation under the assumption 
that the times involved are long compared to the relaxation time r, as denned below. But, 
because the fluctuation-dissipation theorem is really a quantum-mechanical relationship, it 
will tell us something more than just the laws of classical diffusion, which will emerge only 
in the limit h — > O, or equivalently in the limit of long times and high temperatures. In the 
opposite limit where kTAt <C h, the usual linear dependence Ax 2 ~ At, will turn out to 
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give way to a universal behavior Ax 2 ~ In At, which probably should be interpreted as a 
diffusion driven by quantal zero-point motions rather than by thermal kinetic energy. The 
logarithmic behavior will follow from a general formula (12) for (Ax 2 ), which will hold for 
all times long compared to r, given the assumption of constant mobility ji. In what follows 
we will derive this general formula, discuss the limiting cases just alluded to, and show that 
some deviations from classical behavior may be observable on the basis of current experi- 
mental technique. 

In recent years there have been several efforts 1-4 to understand the dynamics of a quantum 
particle coupled to a heat bath. Insofar as our work overlaps those efforts, our results appear 
to agree. The main difference is that the cited papers make far-reaching assumptions about 
the nature of the medium (heat bath) in which the particle moves, and require the coupling 
between particle and bath to be linear (meaning in effect that the coupling is weak). In 
contrast, we only use that the response to a weak external perturbation is linear, allowing 
the coupling of the particle to the bath and/or environment itself to be strong, as it will 
in fact be in most situations. On the other hand we will predict only the mean-square 
displacement, whereas the more special treatments can in principle yield the full density 
operator as a function of At. 



The fluctuation-dissipation theorem, as usually stated, refers to the Fourier transforms of 
the autocorrelation and response functions. Let x(t) be some dynamical variable (operator) 
in the Heisenberg picture, and let f(t) be an infinitely weak external force applied to x 
at time t. (We will not need the more general form of the theorem in which the external 
coupling is to a different variable y.) the response function R(t) is defined by the relation 



THE FLUCTUATION-DISSIPATION 



THEOREM IN THE TIME DOMAIN 




(1) 
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where (•••)/ denotes the expectation value in the presence of the force, assuming the system 
of which a; is a variable to have been in thermal equilibrium with temperature T at early 
times; and (• • -) is the same expectation value for zero force. Also let 

C(t) = ±(x(t)x(0)+x(0)x(t)) 

be the "autocorrelation" or "two-point" function in equilibrium at temperature T. [Or, if you 
prefer, you can subtract off (x(t))(x(0)) = (x(0)) 2 from this definition without invalidating 
what follows. This would be equivalent to working with x — (x) in place of x.] Then the 
fluctuation-dissipation theorems stated in the frequency domain is (with (3 — 1/ kT) 

ImR(v) = h- 1 tanh(7r (3his) ~ C(v). (2) 

[We are using the following definition of fourier transform T = (• • •): 

fa) = J dtrY(t), 

where l x = e 2mx .] 

Our first job is to transform this relation to the time domain. To that end, let us introduce 
in place of R(t) (which vanishes for t < by virtue of causality) the equivalent odd function 

R(t) = sgn(t)R(\t\). 

It is then easy to check that 2iImJ r (_R) = J-'(R), whence (2) can be written in the equivalent 
form: 

?{R) = ^tanh(vr^z/)^(C). (3) 

[In fact it is actually this form, rather than (2), that comes out initially in the most straight- 
forward derivation of the fluctuation-dissipation theorem; it is thus more appropriate to view 
(2) as a consequence of (3) rather than vice versa.] By taking the Fourier transform of (3) 
we could now express R(t) as a convolution of C (t), but our main interest here is to do the 
opposite. Let us therefore solve (3) for C, obtaining 
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C[y) = (-ih/2) coth(irPhv)[F(R)](is) + c5(v), (4) 

where c is a constant and where, for definiteness, the principal part of coth may be taken: 
P(cotha;) = rf/rfxlnsinh \x\. The ambiguity in 1/ tanh(7r/5^z/) is just a term proportional to 
S(u), which would drop out of (4) anyway, since it would be multiplying the odd function 
F{R)] The Fourier transform of (4) reads 

C = l -^F(cothirPkv)*R + c, (5) 

determining C, up to an additive constant, in terms of the Fourier transform 

.F(coth7r/3/M/) = (i/ph) coth^t / (3h) . (6) 

In Eq. (6), the coth on the right-hand side is also to be understood as a principal part, 
but unlike before, this choice is forced on us, because the addition of any 8(t) piece to 
coth(7rt / f3K) would spoil its oddness, in disagreement with the oddness of the left-hand side 
of (6). Understanding all coth's to be principal parts, then, we have finally (in view also of 
the definition of R) the following explicit formula for C(t) in terms of R(t): 

1 r°° 

C(t) = —J_Jt'sgn(t'-t)R(\t f -t\) 

x coth(7rt' / Ph) + c. (7) 

[The appearance of the undetermined constant c is due to the possibility of redefining the 
zero of x without affecting (1). By working with the alternative definition of C(t) mentioned 
just before Eq. (2), we would remove this ambiguity, and correspondingly could set c = 0, 
given some assumptions on the asymptotic behavior of x and R] 

THE MEAN-SQUARE DISPLACEMENT (Ax 2 ) 

Now the mean-square displacement of x due to equilibrium fluctuations in time At is (Ax 2 ), 
where Ax = x(t + At) — x(t). Taking t = for convenience, we have (since the equilibrium 
state is time independent) 
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(Ax 2 ) = ([x(At) - x(0)} 2 ) = (x(At) 2 ) + (x(0) 2 ) 



({x(At),x(0)}) =2C(0) -2C(At), 



or 




(8) 



Combining this result with (7) gives us a general equation for (Ax 2 ) in terms of the response 
function R: 



where for brevity we have set Q = ir//3h. Here, as before, the principal part of the coth is 
to be understood. Notice that the undetermined constant c in (7) has dropped out of this 
result. 



At this stage, let us specialize i to be a Cartesian coordinate of an otherwise free particle 
immersed in a homogeneous medium with temperature T. For an idealized inertialess Brow- 
nian particle, the response to a weak external force would be immediate motion at velocity 
v = ///, \x being the "mobility;" in other words, R would be the step function R(t) = //O(t). 
However this idealization is plainly too unrealistic, because it leads to a divergent result in 
(9). [In this sense we might say that the fluctuation-dissipation theorem knows that parti- 
cles have inertia.] A more reasonable Ansatz for R must incorporate a "relaxation time" or 
"rise time" r representing the time it takes the particle to accommodate itself to any sudden 
change fit). Such an Ansatz is, for example., 



which describes the classical motion of a particle subject to viscous friction. Without making 
so specific a choice, however, we will employ a cruder cutoff which should be adequate for 
times much greater than r: 




cothft(i' -t)}, 



(9) 



QUANTUM BROWNIAN MOTION 



R(t) = /i(l - e-*/ r )0(t) 



(10) 
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R(t) =fj,G(t-T). 



(11) 



With this R, (9) can be integrated exactly [using the distributional identity, P(cothx) = 
d/dxhi sinh \x\ to produce the following fundamental equation of quantum Brownian motion: 



where again Q = ir/(3h. 

Now strictly speaking, there is the inconsistency in our derivation of (12) that C(t) is ill 
defined for a particle moving in an unbounded space, because (x 2 ) in equilibrium would 
be infinite, and (8) would therefore assume the indeterminate form (Ax 2 ) = oo — oo. To 
overcome this problem, one could confine the particle in a very long "box" (confining po- 
tential), it being intuitively clear that this could alter neither (Ax 2 ) nor R(t) in the limit of 
an infinitely large such box. 



The possible limiting cases of (12) are determined by the relative magnitudes of the three 
times r, j3h, and At, which we may call, respectively, the relaxation time, the "quantum 
time," and the "diffusion time." A priori there would be essentially 3!=6 distinct cases, but 
since we must have At ^> r in order to apply (12), we will limit ourselves to only three 
of them. [It is nonetheless instructive to notice that (12) becomes self-contradictory for At 
near r since it then equates an intrinsically positive expression to a negative right-hand 
side. This implies that (11) could not be the exact response function for any system, even 
in principle. More generally, one can derive from (7) and the definition of C(t), a positivity 
criterion which any putative response function must fulfill in order to be physically viable. 
We do not know how restrictive this criterion is in practice, but we have checked that the R 




(12) 



THREE LIMITING CASES OF 



THE GENERAL FORMULA (12) 
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of (10) yields a mean-square displacement which is non-negative for all times, as one might 
have expected.] 

Case 1: < At. This is the classical limit, and (12) reduces to the classical relation 

1 -(Ax 2 ) = ^/P)At = f ,kTAt, (13) 
or fikT(At — r) if the leading correction is retained. 

Case 2: r < At C f3h. This is the extreme quantum limit, in which (time) (energy) <C U 
for the time scale set by the diffusion time At and the energy scale set by the thermal energy 
kT. In this limit (12) reduces to 

W ) = ^, (14) 

Z 71 T 

or 



^ln 




rAr 


2 

- 1 


7T 




. T _ 
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if somewhat more precision is desired. It is noteworthy that the temperature has disap- 
peared entirely from this expression (except insofar as it influences \i and r) suggesting a 
quantum Brownian motion due entirely to "zero-point" fluctuations, which are present even 
at absolute zero. Indeed, the striking logarithmic dependence in (14) could also have been 
derived by first taking the zero-temperature limit of the fluctuation-dissipation theorem it- 
self, and only then applying it to the R of a diffusing particle. 

Case 3: r <C f3h <C At. Intermediate between cases 1 and 2, this situation might be described 
as one in which the relaxation occurs on quantum time-scales, although the diffusion time 
itself is already classically long. [A suggestive way to rewrite the inequality r <^ f3% is as 
the relation between diffusion constants, -D c i assica i <^ -D qU antum, where -D c i assica i = and 
-Dquamtum — h/m, with m taken from the "viscous damping" relation r = jj,m envisaged in 
(10).] In this case (12) reduces to 
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M = ^^^ (15) 

which one can interpret as the result of two-stage spreading which follows the quantum law 
(14) up to the time £q = (3h/2ir, and thereafter continues according to the classical law (13), 
with the second term in (15) remaining forever as a kind of residue of the quantum era. In 
order for this residue to be significant, we need fj,At/P < (^h/n) \n(Ph/2rcT), or 

pn, phi it 
At < — In 

which can occur nontrivially (i.e., without reducing to case 2) only if /3K/t is exponentially 
big, so that \n(ph/r) > At/ph ^> 1. Taking m = r//i as earlier, this amounts to a 
requirement that the particle be extremely light: 

m < {ph/fi)e- wAt/ ^. (16) 

REMARKS AND NUMERICAL ESTIMATES 

Equations (13), (14), and (15) are all special cases of the more general relation (12), which 
should be valid whenever r <C At. In other situations, or for more general response functions 
R(t), one must refer back to (9) itself, from which the spreading can always be computed as 
long as R(t) is known. A particularly interesting response function to treat would be (10), 
and another interesting case might be a particle moving in a superfluid. 



In the zero-temperature limit, i.e., in case 2 above, our formula (14) may be compared with 
a result of Ambegaokar, 4 who used a path-integral formalism, and assumed a linear coupling 
between the particle and an environment comprising an infinite collection of harmonic oscil- 
lators. He obtained an expression for the mean-square displacement of a Brownian particle 
in the quantum regime which corresponds to our result given in (14), if we make certain 
identifications. According to Ambegaokar (with a presumed misprint corrected), 

((Ax 2 )) = [(h/n 2 )/ im ]ln\(t^)\ + const. , (17) 
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for (I/7) < t < ((3%). Here, (Ax 2 ) is given by a density operator p which reduces to a 5 
function at t — 0, and oj c defines an upper frequency cutoff beyond which the linear rela- 
tionship between particle velocity and environmental friction breaks down. Also, judging 
from Eq. (5.3) of Ref. 4, it appears natural to identify I/7 with our r, and therefore l/Wy 
with our /i. If we do so, and also equate oj c to r _1 , then we recover (14) from (17) with the 
constant set to zero. 

In connection with (14) one can ask the following question: Classically, what kind of response 
function would lead to a logarithmic law of diffusion? If we take the H — > limit of (9), we 
find that the relevant response function should be proportional to 1/t, which is physically 
impossible. This implies that the effect described by (14) is of purely quantum- mechanical 
origin. 

Finally, let us estimate the thresholds of time and temperature at which significant devia- 
tions from classical behavior should appear. In order to be in the "pure quantum regime," 
we need At <C (3% which can also be written in the time-energy form, kTAt <C U. Taking 
T ~ 1(T 6 deg (cf Ref 6) and At ~ 1(T 6 sec yields kTAt/h ~ 0.1, which ought to be well 
within the "pure quantum regime," meaning that (14) should apply if the relaxation time 
is short enough (and the "reservoir" in thermal equilibrium). For higher temperatures or 
longer times, deviations of the sort described by (15) might be observable if r is small enough 
and a condition like (16) is satisfied. 
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